Differential directed flow in Au+Au coflisions 
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We present experimental data on directed flow in semi- 
central Au+Au collisions at incident energies from 90 to 
400-A MeV. For the first time for this energy domain, the data 
are presented in a transverse momentum differential way. We 
study the first order Fourier coefficient vi for different par- 
ticle species and establish a gradual change of its patterns 
as a function of incident energy and for different regions in 
rapidity. 

PACS: 25.70.Lm, 21.65.+f, 25.75.Ld 

The main motivation for the study of relativistic heavy 
ion collisions is to learn about the equation of state (EoS) 
of nuclear matter . The collective phenomena |g| (see 
ref. for recent reviews) were proposed as probes 

that preserve information about the transient hot and 
compressed state created in such collisions. The (in- 
plane) directed flow (also called sidewards flow or trans- 
verse flow) was much studied, both experimentally |5|-p^] 
and theoretically ||r^-|2l|. The dependences of this phe- 
nomenon on incident energy, particle type and central- 
ity have been determined experimentally (see ref. [^[^] 
and references therein), particularly for the Au-|-Au sys- 
tem. The balance energy, Ebai, which is the energy of 
"disappearance of flow" , was recently directly measured 
for Au-|-Au collisions It was found to be around 

40-A MeV, somewhat lower than previous extrapolations 
1^,11^ . The main features of the directed flow have been 
reproduced by the theoretical models, but a final conclu- 
sion on the EoS has not yet been achieved. As pointed 
out early on [p"3|-p5|, the momentum dependent interac- 
tions (MDI) play a crucial role in the determination of 
the EoS. Together with the (in- medium) nucleon-nuclcon 
cross section (cr„„), MDI has a marked influence on the 



directed flow |T^JT^,|2l|. Moreover, consistency is needed 
in deriving EoS together with both MDI and (t„„ . 

All the above-mentioned results have been obtained 
in a transverse momentum (pt) integrated way, usually 
by studying the average in-plane transverse momentum, 
(px), as a function of rapidity, y. The analysis of flow 
in terms of pt dependence of the first order Fourier co- 
efficient, which we shall call "differential directed fiow" 
(DDF), has been proposed by Pan and Danielewicz [ p^ , 
who found it sensitive to MDI. Li and Sustich have stud- 
ied the sensitivity of DDF to both EoS and nuclcon- 
nucleon cross section ((T„„) around the balance energy 
They found complex patterns of DDF, a vanish- 
ing slope of the (px) — y distribution being the result 
of averaging over transverse momenta. The DDF was 
extensively exploited at AGS energies both experimen- 
tally, by the E877 collaboration |2^, and theoretically 
psf . It was found to be sensitive to the characteristics 
of the transverse expansion and on the EoS, especially 
at high Pt p^ ]. For our energy domain, apart from its 
importance in revealing the dynamics of the heavy ion 
collisions, the differential flow could ultimately help to- 
wards pinning down the stiffness of the EoS of hot and 
dense nuclear matter formed in such collisions p8| , p^ . 

In this Rapid Communication we present experimen- 
tal data on directed flow for Au-|-Au collisions at incident 
energies of 90, 120, 150, 250 and 400-A MeV. Along with 
the integrated distributions, the data are presented in a 
Pt differential way for the first time for this energy do- 
main. The DDF is analyzed for different particle species 
and as a function of rapidity. 

The data have been measured with a wide phase-space 
coverage using the FOPI detector |^ at GSI Darmstadt. 
The reaction products were identified by charge (Z) in the 
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forward Plastic Wall (PW) at 1.2° < Oiab < 30° using 
time-of-flight and specific energy loss. In the Central 
Drift Chamber (CDC), covering 34° < Oiat < 145°, the 
particle identification is obtained using magnetic rigidity 
and the energy loss. For more details on the detector 
configuration for this experiment see ref. |^^ . 

For the centrality selection we used the charged parti- 
cles multiplicities, classified into five bins. The variable 
Erat = X^i ^-L,i/ Si -^IM (t^i^ sums run over the trans- 
verse and longitudinal cm. kinetic energy components 
of all the products detected in an event) has been ad- 
ditionally used for a better selection of the most central 
collisions. The geometric impact parameters interval for 
our centrality region M4, for which the results are pre- 
sented here, is 2.0-5.3 fm. We note that in this bin the 
directed flow exhibits a maximum as a function of cen- 
trality |,|ll. 

To compare different incident energies, we use normal- 
ized center-of-mass (cm.) transverse momentum (per 
nucleon) and rapidity, defined as: 
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where the subscript P denotes the projectile. 

The reaction plane has been reconstructed event-by- 
event using the transverse momentum method pTf . All 
particles in an event have been used for the reaction plane 
reconstruction, excluding the particle-of-interest to pre- 
vent autocorrelations and a window around midrapidity 
^|y(o)| < 0.3) to improve the resolution. The correction 
of the extracted values due to the reconstructed reaction 
plane fluctuations has been done using the recipe of QUi- 
trault |2^. The correction factors, l/(cosA0), where 
is the resolution of the reaction plane azimuth, are pre- 
sented in Table | for the five incident energies, along with 
the projectile momenta per nucleon in the cm. system, 

Pp I^P- 

TABLE I. Projectile momentum (per nucleon) in the cen- 
ter of mass and the correction factor for the reaction plane 
resolution. 



E (A MeV) 


90 


120 


150 


250 


400 


p'P/Ap (MeV/c) 


205 


237 


265 


342 


433 


l/(cos A(/)) 


1.51 


1.18 


1.10 


1.05 


1.04 



An important ingredient in the present analysis is the 
correction for distortions due to multiple hit losses. As an 
example, in case of PW, despite its good granularity (512 
independent modules p5| ] ) , average multiple hit probabil- 
ities of up to about 9% at 400-A MeV are registered for 
the multiplicity bin M4. Because of the directed flow, the 
average number of particles detected over the full PW 
subdetector is 2 times higher than out of the reaction 
plane. This leads to an underestimation of the directed 
flow and needs to be taken into account. We developed 
a correction procedure based on the experimental data, 
by exploiting the DDF left-right symmetry with respect 



to midrapidity. The correction depends on pf^ and inci- 
dent energy. It reaches 12% for the energy of 400-A MeV 
and is almost negligible at 90-A MeV. The procedure 
was checked and validated using IQMD (Isospin Quan- 
tum Molecular Dynamics |l^) events passed through a 
fuU GEANT simulation of the detector. 
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FIG. 1. Integrated vi values as a function of rapidity for 
five incident energies, M4 centrality bin, for Z=2 particles. 

We characterize the directed flow by the first order 
Fourier coefficient: vi = (cos(</))) where 4> is the angle 
with respect to the reaction plane. In Fig. |^ we present 
for 7j—2 particles the values of the vi coefficient inte- 
grated over transverse momentum. These distributions 
are for the forward hemisphere and were obtained by 
combining the information of both PW and CDC. Note 
that, for the phase space coverage of CDC, Z=2 sam- 
ple comprises only alpha particles. We mention that the 
slopes at midrapidity of the equivalent {px) — y'"' distri- 
butions are in perfect agreement with the FOPI results 
from Phase I data |^ . 

Two important features can be emphasized from the 
distributions presented in Fig. |l|. First, we are in a region 
where the flow is strongly increasing as a function of inci- 
dent energy ||^ , as a result of stronger compression. This 
is contrary to expectations based on ideal hydrodynamics 
p9[ and has been related to a possible liquid-gas phase 
transition which breaks the scaling behaviour |^,^ . Sec- 
ond, the shape as a function of rapidity is evolving with 
the incident energy. For higher energies, a saturation 
of the vi coefficient is seen for the projectile rapidities. 
Weighted with the (rapidity dependent) transverse mo- 
mentum spectra it leads to the well known so-called "S 
curve" in the {px) — y distributions P,|7|-|^,pT[| which was 
interpreted as a result of "bounce-off" from the specta- 
tors. At lower energies, the steady evolution of vi as a 
function of y points to a poorer separation of the partic- 
ipant and spectator contributions, probably as a result 
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of a relatively stronger Fermi motion. Intuitively, longer 
passing times might as well allow for larger cross-talk and 
"equilibration" between participant and spectators. 
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FIG. 2. The dependence of the vi coefficient on transverse 
momentum for 250 -A MeV incident energy, M4 centrality bin, 
for different particles. The lines are polynomial fits to guide 
the eye. The arrows mark the values of the average 'for 
each particle. 

In Fig. H we show the differential flow for parti- 
cles with charge Z=l to Z=6 for the incident energy 
E=250-A MeV. We selected the rapidity window 0.5< 
y{o) <o.7, corresponding to the maximum of the directed 
flow in the {px)-y representation. From here on, to keep 
a clean particle selection on Z, we limit our pt domain 
to the PW acceptance (which translates in a rapidity- 
dependent upper pf himit). The error bars on vi include 
both statistical and systematical errors, the latter com- 
ing from the apparatus acceptance. The horizontal bars 
indicate the pf'^hin width. The arrows, labeled by the 
corresponding symbols, are indicating the average values 
of P(''''for the fragments analyzed. The lines here and in 
the following figures are polynomial fits to guide the eye. 
For the chosen rapidity window the multiplicities of the 
fragments up to Z=6 relative to Z=l are in the ratios 
1000/380/68/20/14/8. We first remark that the flow of 
heavier fragments is larger, a fact established in an inte- 
grated way by many experiments and interpreted 
as a signature of a collective phenomenon |^^. As was 
demonstrated in ref. this is the result of the inter- 
play of the thermal and collective motion. For the lighter 
particles the thermal contribution is higher relative to the 
collective one, leading to a lower "apparent" flow. One 
can see in Fig. || that not only the magnitude is increas- 
ing, but also that the pattern of vi changes for heav- 
ier particles, for which a maximum develops at steadily 
lower scaled momenta. The increase of the anisotropy as 
a function of pt is very rapid and compatible with a linear 



one only in a very narrow pt region. Note that at AGS 
energies a linear dependence of vi on pt was established 
experimentally for different particle species over a broad 
Pt interval P2[ . 
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FIG. 3. The dependence of the vi coefficient on transverse 
momentum for five incident energies, M4 centrality bin, for 
Z=l and Z=2 particles. 

Fig. H presents, for our five incident energies, the de- 
pendence of the vi coefficient on the normalized trans- 
verse momentum for Z=l and Z=2 particles in the ra- 
pidity window 0.5< y'^^^ <0.7. The two particle species 
show different patterns of DDF as function of incident 
energy. For the lowest incident energy, 90-A MeV, while 
for Z=l particles the flow is close to zero for the low 
pf ^part and shows a later onset, in the case of Z=2 frag- 
ments one can notice a steady increase of vi as a func- 
tion of pf\ Although with smaller magnitude (as we 
are well above Ei,ai), the Z=l case is similar to the be- 
havior predicted by Li and Sustich at Ebai ||2|]. The 
coexistence of two flow patterns, attractive and repul- 
sive, was explained in ref. ||2l[ by the fact that the low 
Pt particles suffer more the influence of the attractive 
mean field, leading to a negative vi, while at high pt the 
repulsive nucleon-nucleon scatterings are responsible for 
the positive flow. The particular feature of DDF seen at 
90-A MeV is rapidly changing as a function of incident 
energy, the increase of vi as a function of pt becoming 
gradually steeper for higher energies. 

For both particle species there is a gradual develop- 
ment of a limiting value of DDF at high pf^as the in- 
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cident energy increases. Part of these high-pt parti- 
cles could have been emitted at a pre-equilibrium stage, 
therefore not reaching the maximum compression stage 
of the reaction. For the corresponding rapidity window, 
the average transverse momentum of Z=:l particles has 
a more pronounced variation compared to Tj—1 (see Ta- 
ble U), presumably as a result of a different participa- 
tion of different particle types in the reaction dynamics 
at these energies Also, Z=2 particles can originate 
to a large extent from secondary decays of heavier frag- 
ments [|l9|. It would be important to study these details 
using transport models. 



TABLE II. Average normalized transverse momentum 
for particles with Z=l and Z=2 for the rapidity window 
=0.5-0.7. 



Energy (A MeV) 


90 


120 


150 


250 


400 


(pD z=i 


0.74 


0.71 


0.67 


0.63 


0.60 




0.51 


0.50 


0.49 


0.48 


0.48 
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FIG. 4. vi as a function of transverse momentum in three 
windows in rapidity for particles with Z=l and Z=2 for the 
incident energy of 90-A MeV. 

From our differential flow values at 90-A MeV one can 
infer that Etai is not a single-defined value but is de- 
pendent on particle type and on transverse momentum. 
We have observed similar trends in what concerns Etran, 
the energy of the transition from in-planc to out-of-plane 
azimuthal enhancement at midrapidity p6| . As Ebai is 



decreasing as a function of pt, its integrated value can- 
not be compared among various experiments unless their 
different thresholds are taken into account. Moreover, 
Ehai can depend on rapidity as well, as shown in Fig. ^, 
where the DDF is presented for Z—1 and Z—2 particles 
for three windows in rapidity for the incident energy of 
90-A MeV. Notice again the difference between the two 
particle species in what concerns the patterns of the flow 
and that at projectile rapidities (triangles) the anisotropy 
of Z=2 particles already develops a maximum, in contrast 
to Z=l particles. 
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FIG. 5. Same as Fig. H, but for the incident energy of 
400-A MeV. 

While at our lowest incident energy of 90-A MeV the 
DDF may help towards pinning down the in-medium 
nucleon-nucleon cross section pl| , at higher energies it 
may bear significance on the MDI and ultimately on 
the EoS. The rapidity dependence of the DDF at the in- 
cident energy of 400-A MeV is presented in Fig. ||for Z=l 
and Z=2 particles. Again, at all rapidities there is a clear 
difference in shape between the two particle species. For 
Z=2 particles notice that at projectile rapidities (y^'^^=l) 
the anisotropy is highest (see also Fig. while the (px ) 
values would already show a decreasing trend . 

For both incident energies. Fig. I and |, one can see 
that the shape of the distributions varies as a function 
of rapidity, giving additional support for the influence of 
the dynamics of the collision on the flow values. Also, the 
maximum of vi for Z=2 (for ?/*-'^^=0. 9-1.1) for the energy 
of 400-A MeV occurs at a different value of transverse 
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momentum compared to 90-A MeV, either in terms of 
scaled or not scaled momenta. 

In summary, we have shown first experimental results 
on differential directed flow in heavy ion collisions at in- 
cident energies from 90 to 400-A MeV. Our pt differen- 
tial study allows us to show that not only is the fiow 
larger, but also its pattern is different for heavier parti- 
cles. In addition, it has a complex evolution as a function 
of beam energy. The differential flow pattern of Z=l and 
Z=2 particles was studied as well for different regions in 
rapidity. 

It is expected that the DDF will help to clear 
the difficulties (and sometimes inconsistencies) noticed 
in previous studies devoted to directed flow using 
comparisons with transport models, either of QMD- 
type [^|6|jl|J|,|l^ or from the Boltzmann-Uehling- 
Uhlenbeck (BUU) family [|l| ,^8|,^o| . Present comparisons 
of the experimental differential directed ffow to BUU ||3|] 
or IQMD models show sensitivities to the EoS. But, 
if the EoS is to be extracted from such types of compar- 
isons, more effort is needed, especially in disentangling 
the contributions of momentum dependent interactions 
and in-medium nucleon-nucleon cross sections from the 
EoS itself. It is in the energy range of the present study 
that both MDI and <Jnn have a very important contri- 
bution and cannot at all be neglected. Along with the 
centrality dependence of directed flow [ p"3HT^JlO[ | and with 
the differential elliptic flow ||2^ , the DDF can contribute 
to a good extent in fixing these quantities . Prop- 

erly incorporating the fragment production [|9[ and a 
consistent treatment of the momentum dependent inter- 
actions 1 17 19| , p0[| and in-medium nucleon-nucleon cross 
sections arc essential prerequisites for this task. 
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